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We develop a field theoretical approach to the cold interstellar medium (ISM). We show that a 
non-relativistic self-gravitating gas in thermal equilibrium with variable number of atoms or frag- 
ments is exactly equivalent to a field theory of a single scalar field 4>(x) with exponential self- 
interaction. We analyze this field theory perturbatively and non-perturbatively through the renor- 
malization group approach. We show scaling behaviour (critical) for a continuous range of the 
temperature and of the other physical parameters. We derive in this framework the scaling relation 
AM (R) ~ R dH for the mass on a region of size R, and Av ~ R q for the velocity dispersion where 
q = ^ (riff — 1). For the density-density correlations we find a power-law behaviour for large dis- 
tances ~ |fi — r*2| 2dff_6 . The fractal dimension dn turns to be related with the critical exponent 
v of the correlation length by dn = 1/v. The renormalization group approach for a single com- 
ponent scalar field in three dimensions states that the long-distance critical behaviour is governed 
by the (non-perturbative) Ising fixed point. The corresponding values of the scaling exponents are 
v — 0.631..., dn = 1.585... and q = 0.293.... Mean field theory yields for the scaling exponents 
v — 1/2, dn = 2 and q = 1/2. Both the Ising and the mean field values are compatible with the 
present ISM observational data: 1.4 < du < 2, 0.3 < q < 0.6 . 

As typical in critical phenomena, the scaling behaviour and critical exponents of the ISM can be 
obtained without dwelling into the dynamical (time-dependent) behaviour. 

The relevant role of selfgravity is stressed by the authors in a Letter to Nature, September 5, 1996. 
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I. INTRODUCTION AND RESULTS 



The interstellar medium (ISM) is a gas essentially formed by atomic (HI) and molecular (H2) hydrogen, distributed 
in cold (T ~ 5 — 50K) clouds, in a very inhomogeneous and fragmented structure. These clouds are confined in 
the galactic plane and in particular along the spiral arms. They are distributed in a hierarchy of structures, of 
observed masses from 1 Mq to 10 6 M Q . The morphology and kinematics of these structures are traced by radio 
astronomical observations of the HI hyperfine line at the wavelength of 21cm, and of the rotational lines of the CO 
molecule (the fundamental line being at 2.6mm in wavelength), and many other less abundant molecules. Structures 
have been measured directly in emission from O.Olpc to lOOpc, and there is some evidence in VLBI (very long based 
interferometry) HI absorption of structures as low as 10~ 4 pc = 20 AU (3 10 14 cm). The mean density of structures is 
roughly inversely proportional to their sizes, and vary between 10 and 10 5 atoms / 'cm 3 (significantly above the mean 
density of the ISM which is about 0.1 atoms/cm 3 or 1.6 10 -25 g/cm 3, ). Observations of the ISM revealed remarkable 
relations between the mass, the radius and velocity dispersion of the various regions, as first noticed by Larson [0, 
and since then confirmed by many other independent observations (see for example ref. Q). From a compilation of 
well established samples of data for many different types of molecular clouds of maximum linear dimension (size) R, 
mass fluctuation AM and internal velocity dispersion Av in each region: 

AM(R) ~ R dH , Av ~ R q , (1) 

over a large range of cloud sizes, with 10 -4 — 10~ 2 pc < R < 100 pc, 

1.4 < d H < 2, 0.3 < q < 0.6 . (2) 

These scaling relations indicate a hierarchical structure for the molecular clouds which is independent of the scale 
over the above cited range; above 100 pc in size, corresponding to giant molecular clouds, larger structures will be 
destroyed by galactic shear. 
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These relations appear to be universal, the exponents g?#, q are almost constant over all scales of the Galaxy, and 
whatever be the observed molecule or element. These properties of interstellar cold gas are supported first at all from 
observations (and for many different tracers of cloud structures: dark globules using 13 CO, since the more abundant 
isotopic species 12 CO is highly optically thick, dark cloud cores using HCN or CS as density tracers, giant molecular 
clouds using 12 CO, HI to trace more diffuse gas, and even cold dust emission in the far-infrared). Nearby molecular 
clouds are observed to be fragmented and self-similar in projection over a range of scales and densities of at least 10 4 , 
and perhaps up to 10 6 . 

The physical origin as well as the interpretation of the scaling relations ([!]) are not theoretically understood. The 
theoretical derivation of these relations has been the subject of many proposals and controversial discussions. It is 
not our aim here to account for all the proposed models of the ISM and we refer the reader to refs. Q for a review. 

The physics of the ISM is complex, especially when we consider the violent perturbations brought by star formation. 
Energy is then poured into the ISM either mechanically through supernovae explosions, stellar winds, bipolar gas flows, 
etc.. or radiatively through star light, heating or ionising the medium, directly or through heated dust. Relative 
velocities between the various fragments of the ISM exceed their internal thermal speeds, shock fronts develop and 
are highly dissipative; radiative cooling is very efficient, so that globally the ISM might be considered isothermal 
on large-scales. Whatever the diversity of the processes, the universality of the scaling relations suggests a common 
mechanism underlying the physics. We propose that self-gravity is the main force at the origin of the structures, 
that can be perturbed locally by heating sources. Observations are compatible with virialised structures at all scales. 
Moreover, it has been suggested that the molecular clouds ensemble is in isothermal equilibrium with the cosmic 
background radiation at T ~ 3K in the outer parts of galaxies, devoid of any star and heating sources Q. This colder 
isothermal medium might represent the ideal frame to understand the role of self-gravity in shaping the hierarchical 
structures. Our aim is to show that the scaling laws obtained are then quite stable to perturbations. 

Till now, no theoretical derivation of the scaling laws eq. (Q) has been provided in which the values of the exponents 
are obtained from the theory (and not just taken from outside or as a starting input or hypothesis). 

The aim of these authors is to develop a theory of the cold ISM. A first step in this goal is to provide a theoretical 
derivation of the scaling laws eq.(Jl]), in which the values of the exponents du , q are obtained from the theory. For 
this purpose, we will implement for the ISM the powerful tool of field theory and the Wilson's approach to critical 
phenomena |f2| . 

We consider a gas of non-relativistic atoms interacting with each other through Newtonian gravity and which are 
in thermal equilibrium at temperature T. We work in the grand canonical ensemble, allowing for a variable number 
of particles N. 

Then, we show that this system is exactly equivalent to a field theory of a single scalar field 4>(x) with exponential 
interaction. We express the grand canonical partition function Z as 



V<p e 



s[<t>(-)] 



(3) 



where 



s[4>(.)] 



T, 



d A x 



T eff = 4tt 



eff 

Gm 2 



T 



-(V0) 2 -p 2 e^ 



V = 



Gm 7/2 VT , 



(4) 



m stands for the mass of the atoms and z for the fugacity. We show that in the <^>-field language, the particle density 
expresses as 



< p{r) >- 



I 



< V>(r) >= 



< e"' 1 > 



(5) 



T eff T eff 

where < . . . > means functional average over <p{ ) with statistical weight e~ s ^'". Density correlators are written as 
C(f 1 ,f 2 ) = < p{fx)p{r 2 ) > - < pin) >< p(f 2 ) > 



The 0-field defined by eqs. 



T, 



eff 



< e^ Fl) > - < e^ (?l) > < > 



(6) 



has remarkable properties under scale transformations 
x —>■ x\ = Xx , 
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where A is an arbitrary real number. For any solution 4>{x) of the stationary point equations, 

V 2 0(f) + /i 2 e*™ = , (7) 
there is a family of dilated solutions of the same equation (j?]), given by 

<Px(x) = (j)(\x) +logA 2 . 

In addition, S[4>\(.)] = \ 2 ~ D S[</>(.)]. 

We study the field theory (||)-(|4|) both perturbatively and non-perturbatively. 

The computation of the thermal fluctuations through the evaluation of the functional integral eq. ([|) is quite non- 
trivial. We use the scaling property as a guiding principle. In order to built a perturbation theory in the dimensionless 
coupling g = ^//iT e jj we look for stationary points of eq.(Q). We compute the density correlator eq.(^|) to leading 
order in g. For large distances it behaves as 

m,K>»-V~ a <.*_ w +0(K-V*) . (8) 

We analyze further this theory with the renormalization group approach. Such non-perturbative approach is the 
more powerful framework to derive scaling behaviours in field theory |l2|~|l4|] . 
We show that the mass contained in a region of volume V = R 3 scales as 

r R K 

<M(R)>=m / < e^ (5) > d 3 x~ mVa + m R» + . . . , 

J \-a 

and the mass fluctuation, (AM(i?)) 2 =< M 2 > — < M > 2 , scales as 

AM (R) ~ R dH . 

Here v is the correlation length critical exponent for the ^-theory (||) and a and K are constants. Moreover, 

< p(r) >= ma + m j- r r~~ 3 for r of order ~ R . (9) 

iirvyl — a) 

The scaling exponent v can be identified with the inverse Haussdorf (fractal) dimension dn of the system 

d H = - . 

v 

In this way, AM ~ R d " according to the usual definition of fractal dimensions . 

From the renormalization group analysis, the density-density correlators (|^) result to be, 

<7(fi,f 2 )~ \n~r 2 \i- e . (10) 

Computing the average gravitational potential energy and using the virial theorem yields for the velocity dispersion, 

Av-R^- 1 ^ . 

This gives a new scaling relation between the exponents djj and q 

q= U--A =l(d H -l). 
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The perturbative calculation (B) yields the mean field value for v |n|. That is, 

v = i , d H = 2 and q = i . (11) 

We find scaling behaviour in the ^-theory for a continuum set of values of fj? and T e *f. The renormalization 
group transformation amounts to replace the parameters /i 2 and T e ff in /3 H and 5 I [0(.)] by the effective ones at the 
scale L in question. 
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The renormalization group approach applied to a single component scalar field in three space dimensions indicates 
that the long distance critical behaviour is governed by the (non-perturbative) Ising fixed point fl2]-fllf| . Very probably, 
there are no further fixed points (It]] . The scaling exponents associated to the Ising fixed point are 

v = 0.631... , d H = 1.585... and q = 0.293... . (12) 

Both the mean field (|ll|) and the Ising (|l|) numerical values are compatible with the present observational values 
(I)-®- 

The theory presented here also predicts a power-law behaviour for the two-points ISM density correlation function 
(see eq.(pX|), 2du — 6 = —2.830 . . ., for the Ising fixed point and 2du — 6 = —2 for the mean field exponents), that should 
be compared with observations. Previous attempts to derive correlation functions from observations were not entirely 
conclusive, because of lack of dynamical range (^3|, but much more extended maps of the ISM could be available 
soon to test our theory. In addition, we predict an independent exponent for the gravitational potential correlations 
(~ where T]i S i ng — 0.037. . . and J] mean field = p3|), which could be checked through gravitational lenses 

observations in front of quasars. 

The mass parameter \i [see eq.(Q)] in the ^-theory turns to coincide at the tree level with the inverse of the Jeans 
length 

1 

~d~j ■ 

We find that in the scaling domain the Jeans distance dj grows as < dj >^ R. This shows that the Jeans distance 
scales with the size of the system and therefore the instability is present for all sizes R. Had dj being of order larger 
than R, the Jeans instability would be absent. 

The gravitational gas in thermal equilibrium explains quantitatively the observed scaling laws in the ISM. This fact 
does not exclude turbulent phenomena in the ISM. Fluid flows (including turbulent re gimes) are probably relevant in 
the dynamics (time dependent processes) of the ISM. As usual in critical phenomena |l2j,|l3[, the equilibrium scaling 
laws can be understood for the ISM without dwelling with the dynamics. A further step in the study of the ISM will 
be to include the dynamical (time dependent) description within the field theory approach presented in this paper. 

If the ISM is considered as a flow, the Reynolds number ReisM on scales L ~ lOOpc has a very high value of the 
order of 10 6 . This led to the suggestion that the ISM (and the universe in general) could be modelled as a turbulent 
flow H . (Larson []lj first observed that the exponent in the power- law relation for the velocity dispersion is not greatly 
different from the Kolmogorov value 1/3 for subsonic turbulence). 

It must be noticed that the turbulence hypothesis for the ISM is based on the comparison of the ISM with the 
results known for incompressible flows. However, the physical conditions in the ISM are very different from those 
of incompressible flows in the laboratory. (And the study of ISM turbulence needs more complete and enlarged 
investigation than those performed until now based in the concepts of flow turbulence in the laboratory). Besides the 
facts that the ISM exhibits large density fluctuations on all scales, and the observed fluctuations are highly supersonic, 
(thus the ISM can not viewed as an 'incompressible' and 'subsonic' flow), and besides other differences, an essential 
feature to point out is that the long-range self-gravity interaction present in the ISM is completely absent in the 
studies of flow turbulence. In any case, in a satisfactory theory of the ISM, it should be possible to extract the 
behaviours of the ISM (be turbulent or whatever) from the theory as a result, instead to be introduced as a starting 
input or hypothesis. 

This paper is organized as follows. In section II we develop the field theory approach to the gravitational gas. A 
short distance cutoff is naturally present here and prevents zero distance gravitational collapse singularities (which 
would be unphysical in the present case). Here, the cutoff theory is physically meaningful. The gravitational gas is 
also treated in a D-dimensional space. 

In section III we study the scaling behaviour and thermal fluctuations both in perturbation theory and non- 
perturbatively (renormalization group approach) . g = fx T e f f acts as the dimensionless coupling constant for the 
non-linear fluctuations of the field (j>. We show that these fluctuations are massless and that the theory scales (behaves 
critically) for a continuous range of values n 2 T e ff. Thus, changing fi 2 and T e // keeps the theory at criticality. The 
renormalization group analysis made in section III confirm such results. We also treat (sect. III.E) the two dimensional 
case making contact with random surfaces and their fractal dimensions. 

Discussion and remarks are presented in section IV. External gravity forces to the gas like stars are shown not to 
affect the scaling behaviour of the gas. That is, the scaling exponents q, du are solely governed by fixed points and 
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hence, they are stable under gravitational perturbations. In addition, we generalize the ^-theory to a gas formed by 
several types of atoms with different masses and fugacities. Again, the scaling exponents are shown to be identical to 
the gravitational gas formed of identical atoms. 

The differences between the critical behaviour of the gravitational gas and those in spin models (and other statistical 
models in the same universality class) are also pointed out in sec. IV. 



II. FIELD THEORY APPROACH TO THE GRAVITATIONAL GAS 



Let us consider a gas of non-relativistic atoms with mass m interacting only through Newtonian gravity and which 
are in thermal equilibrium at temperature T = /3 _1 . We shall work in the grand canonical ensemble, allowing for a 
variable number of particles N. 

The grand partition function of the system can be written as 



00 _JV 

z= y — 



N=0 



N 



Y[ d3pi ^ e~? HN 



1=1 



(2tt) 3 



where 



N 2 1 



7 = 1 



1<7^7<7V 



G is Newton's constant and z is the fugacity. 

The integrals over the momenta pi, (1 < I < N) can be performed explicitly in eq.(Q3) using 



d 3 p 0P 2 
e 2 ™ = 



(27T) 



m 
2tt/3 



3/2 



We thus find, 



(13) 



(14) 



OO 1 

z = y - 



N=0 



m 

2^/3 



3/2 



N 



J 7 = 1 



We proceed now to recast this many-body problem into a field theoretical form P,^, OJlg| . 
Let us define the density 



N 



3=1 



such that, we can rewrite the potential energy in eq.(p_5|) as 



\f3Gm 2 V - 1 - , = l(3Gm 2 [ 

2 i<t^< N \ qi ~ q i\ 2 J\3-iJ\>a 



d 3 x d 3 y 
\x- y\ 



p(x)p{y) . 



(15) 



(16) 



(17) 



The cutoff a in the r.h.s. is introduced in order to avoid self-interacting divergent terms. However, such divergent 
terms would contribute to Z by an infinite multiplicative factor that can be factored out. 
By using 



1 



f - y\ 



= — 47r S(x — y) 



and partial integration we can now represent the exponent of the potential energy eq. ( [lt| ) as a functional integral 0] 



V£ e 



-§ J d 3 x (V^) 2 + Im^itGP J d 3 x p(x) 



(18) 



Inserting this expression into eq.(|l5[) and using cq.(|16[) yields 
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AH 
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i=i 



^ AT! 



W/3J 



N 
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2mJitGI3 £(q) 



N 



X>f e 



•/ d 3 x i(v?) 2 



3/2 2 mv /^G? {(sf) 



It is convenient to introduce the dimensionless field 



Then, 



where 



Z = 



(j){x) = 2m v / 7rG/3 . 



7T J 



Gm 2 



(19) 



(20) 



(21) 



(22) 



The partition function for the gas of particles in gravitational interaction has been transformed into the partition 
function for a single scalar field 4>(x) with local action 



s[<K.)] 



i 



T, 



eff 



d 3 3 



^(V0) 2 ~H 2 e^ 



(23) 



The field exhibits an exponential self-interaction — /i 2 e^ x \ 

Notice that the effective temperature T e f / for the 0- field partition function turns out to be inversely proportional 
to T whereas the characteristic length behaves as ~ T -1 / 4 . This is a duality-type mapping between the two 
models. 

It must be noticed that the term — fj 2 makes the </>- field energy density unbounded from below. Actually, 

the initial Hamiltonian (|l^) is also unbounded from below. This unboundness physically originates in the attractive 
character of the gravitational force. Including a short-distance cutoff [see sec. 2A, below] eliminates the zero distance 
singularity and hence the possibility of zero-distance collapse which is unphysical in the present context. We therefore 
expect meaningful physical results in the cutoff theory Moreover, assuming zero boundary conditions for <p(r) at 
v — ► oo shows that the derivatives of <j> must also be large if is large. Hence, the term i(V(/>) 2 may stabilize the 
energy. 

The action (^) defines a non-renormalizable field theory for any number of dimensions D > 2 [see eq.(^)]. This 
is a further reason to keep the short-distance cutoff non-zero. 
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Let us compute now the statistical average value of the density p(f) which in the grand canonical ensemble is given 



< p(r) > 



OO , 



N=0 



m 

2^/3 



3/2' 



N 



i=i 



(24) 



As usual in the functional integral calculations, it is convenient to introduce sources in the partition function ( pl| ) 
in order to compute average values of fields 



Z[J(.)\ z 

The average value of <f)(f) then writes as 



< <f>(r) >- 



SlogZ 
SJ{r) 



(25) 



(26) 



G 



In order to compute < p(f) > it is useful to introduce 



V[J(.)] = \f3Gm 2 / + J(x) ] [p(y) + J(y) 

2 J\S~y\>a \X-y\ 



(27) 



Then, we have 



A „V[0] 



p(r) e 



Teff r \6J(r) 
By following the same steps as in eqs.(|l8|)-([l9|), we find 



.1=0 



< P(r) > 



1 



Teff r \ SJ(r) ^ N\ 
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/ \ 3/2' 



JY 



JV 



p£ e -/d 3 x[i (VC) 2 -2m v ^3e(x) J(S)] /" / 'jTd 3 g, gfrW^GP E^x^ 

J •* 1=1 

= -T~V 2 J J ^- logZ[J(.)])\ J=0 . 
T eff \SJ[r) J 

Performing the derivatives in the last formula yields 

1 1 r v<f> v 2 ^ e "^7 J**[i<y*r-s* 



|J=0 



< P(r) >- 



T, 



eff 



(28) 



(29) 



One can analogously prove that p(r) inserted in any correlator becomes ~t~J V 2 */^) m the 0- field language. 
Therefore, we can express the particle density operator as 

1 



p(r) 



'eff 



V 2 0(f) 



(30) 



Let us now derive the field theoretical equations of motion. Since the functional integral of a total functional 
derivative identically vanishes, we can write 



V<\> 



SS 



J(r) 



-S[d>(.)]+ J d 3 x J(S) <f,(g) _ Q 



We get from eq.(E3| 



Thus, setting J(r) ee 0, 



V 2 </.(r) + p 2 e*W 



SS _ 1 
W)~~T7f~ f 

< V 2 cj)(f) > + p 2 < >= 



Now, combining eqs.(BOJ) and (|3f|) yields 



< p(r) >= J— < e*& > 



eff 



(31) 



(32) 



By using eq.(pG), the gravitational potential at the point 



U(r) = -Gm 



\x — r\ 



can be expressed as 
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U(r) = -— <t>(r) . 
m 

We can analogously express the correlation functions as 

C(fi,f 2 ) = < p(ri)p(r 2 ) > - < p(n) >< p(r 2 ) > 

2 



(33) 



1 



eff 



|J=0 • 



This can be also written as 



T, 



eff 



< e 



> - < e^' 1 ' > < > 



(34) 



(35) 



A. Short distances cutoff 



A simple short distance regularization of the Newtonian force for the two-body potential is 

Gm 2 



v a {r) 



[l-6(a-r)] , 



9{x) being the step function. The cutoff a can be chosen of the order of atomic distances but its actual value is 
unessential. 

The iV-particle regularized Hamiltonian takes then the form 



N v 2 1 

1=1 1<1,3<N 



(36) 



Notice that now we can include in the sum terms with I = j since v a (0) = 0. 

The steps from eq.(|l4|) to eq.(pl|) can be just repeated by using now the regularized v a (r). Notice that we must use 
now the inverse operator of v a (r) instead of that of 1/r, [— |^V 2 ] , previously used. 

We now find, 



i. e. 



Z a = 



Sa[<K-)] = 



T, 



eff 



d 3 x 



-4>K a 4> - p 2 e+& 



(37) 



(38) 



where K a is the inverse operator of v a , 



K a cj)(r) = I K a (f- f ') 4>{r') d 3 r' 



K a (r-r") £ 
K a (r) admits the Fourier representation, 



d' r" = 8{r — r') 



K a {r) = V.P. 



d 3 p p 2 e l ' y - p 



(27r) 3 cospo 



Actually, K a {f) = for r ^ 0. K a (r) has the following asymptotic expansion in powers of the cutoff a 2 



and then 



K a {r) = -V 2 <5(f) + - V 4 <Kf) + 0(a 4 ) , 



S a y>(.)] = S[4>{.)} + y j d 3 x (V 2 0) 2 + 0(a 4 ) 



(39) 



(40) 



As we see, the high orders in a 2 are irrelevant operators which do not affect the scaling behaviour, as is well known 
from renormalization group arguments. For a — > 0, the action is recovered. 
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B. D-dimensional generalization 

This approach generalizes immediately to .D-dimensional space where the Hamiltonian ( |l4| ) takes then the form 

N 2 i 

H ^Y.h- Gm2 E \a,-a\^ i0T °^ 2 (41) 

1=1 Zm l<l<j<N lQl W 



and 



N 2 1 

^=ES:- Gm2 E ^u-^' at ^ 2 ' ( 42 ) 



2m ' la) — 0,1 

1 = 1 l<l<j<N ,Hl 

The steps from eq.((l|) to (|l]) can be trivially generalized with the help of the relation 

1 



in £)-dimensions and 



at D = 2. 
Here, 



- ^n-2=-C D 8{x-y) (43) 
x — y\ u z 



V 2 \og-J— = -C 2 5{x-y) 
\x-y\ 



2tt°/ 2 

r 

We finally obtain as a generalization of eq.(EJ 



C D = (D - 2) — w - for D ^ 2 and C 2 = 2tt . (44) 



2 -p 2 e*< £ )l 



Z = J J V<j>e T *ff J " * J , (45) 

where 

A^^^^m^/ 2 ^/- 1 , T e// =C C ^. (46) 

We have then transformed the partition function for the Z?-dimensional gas of particles in gravitational interaction 
into the partition function for a scalar field (f> with exponential interaction. The effective temperature T e ff for the 
</!>-ficld partition function is inversely proportional to T for any space dimension. The characteristic length /i -1 
behaves as - T-^- 2 )/ 4 . 

III. SCALING BEHAVIOUR 

We derive here the scaling behaviour of the <j) field following the general renormalization group arguments in the 
theory of critical phenomena |T^ , |l3| ] 

A. Classical Scale Invariance 

Let us investigate how the action ( p3| ) transforms under scale transformations 

x — > x\ = Xx , (47) 

where A is an arbitrary real number. 

In D-dimensions the action takes the form 
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1 eff 



i(V0) 2 - M 2 e^) 



(48) 



Wc define the scale transformed field 4>\(x) as follows 

A (f) = 0(Af)+logA 2 . (49) 

Hence, 

(V^(a?)) 2 = A 2 (V^^a)) 2 , e **(2) = A 2 e*^ 
We find upon changing the integration variable in eq. (^) from x to x\ 

S[M-)] = A 2 " D S^(.)] (50) 

We thus see that the action p^) scales under dilatations in spite of the fact that it contains the dimensionful 
parameter /x 2 . This remarkable scaling property is of course a consequence of the scale behaviour of the gravitational 
interaction (|4l|). 

In particular, in D = 2 the action ( f48|) is scale invariant. In such special case, it is moreover conformal invariant. 

The (Noether) current associated to the scale transformations (|47]) is 

Ji(x) = Xj Tij(x) + 2 Vi<j>(x) , (51) 

where Tij(x) is the stress tensor 

Tij(x) = Vi4>{x) Vj(/>(x) - 5ij L 
and L = i(V0) 2 — /i 2 stands for the action density. That is, 



Ji(x) = (x.V<t> + 2) V l (f>(x) - x, 



By using the classical equation of motion (|52|), we then find 

ViJi{x) = (2-D)L . 

This non-zero divergence is due to the variation of the action under dilatations [eq. (|5C 
If <j>(x) is a stationary point of the action ([l8|) : 

V 2 0(x) + m 2 e 0(S) =0, (52) 
then (j>\(x) [defined by eq.([49|)] is also a stationary point: 

V 2 A (^)+M 2 e**® = Q. 
A rotationally invariant stationary point is given by 

c (r) = logfc^. (53) 



This singular solution is invariant under the scale transformations (|49|). That is 

4>l(r) = 4> c {r) • 

Eq.([53"l) is dilatation and rotation invariant. It provides the most symmetric stationary point of the action. Notice 
that there are no constant stationary solutions besides the singular solution 0o = — oo. 

The introduction of the short distance cutoff a, eq. (|36|) , spoils the scale behaviour (|50|). For the cutoff theory from 
eqs.(||) and (|7|)-([49]), we have instead 

S a [<f> x {.)]=X 2 - D S Xa {<f>{.)] ■ 

For r ~ a, eq.(|53]) does not hold anymore for the spherically symmetric solution (jf{r). For small r and a, using 
eqs.(§§-(i0|) we have 

c M'=°-^ + O(rW). (54) 
That is, 4> c (r) is regular at r = in the presence of the cutoff a. 
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B. Thermal Fluctuations 



In this section we compute the partition function eqs . (J2 1|) and (25) by saddle point methods. 
Eq.(p2|) admits only one constant stationary solution 



4>q = — oo . (55) 

In order to make such solution finite we now introduce a regularization term e fi 2 <p(x) with e < < 1 in the action 
S [eq.(pq)]. This corresponds to an action density 



L = i(V0) 2 



(56) 



where 



u{4>) = —p? e* {s) + e p 2 <j>{x) 



This extra term can be obtained by adding a small constant term — e fj?/T e ff to p(x) in eqs.(|l6|) - (ft8|). This is a 
simple way to make 4>q finite. 

We get in this way a constant stationary point at cf>o = loge where u'((f>o) = 0. However, scale invariance is broken 
since u"(4>o) — — e p 2 ^ 0. We can add a second regularization term to ^8 p 2 <fi(x) 2 to L, (with 5 << 1) in order to 
enforce u"((j)o) = 0. This quadratic term amounts to a long-range shielding of the gravitational force. We finally set: 



u(4>) = -n 2 



MS) 



1 



< ■"' - e <j>{x) - - S (j){xY 



where the two regularization parameters e and 5 are related by 

e(S)=6[l-\o g 5} 

and the stationary point has the value 



Expanding around 0o 



log 5 ■ 



where g = W ji D 2 T e f f and x(x) is the fluctuation field, yields 



— L = o V ^) 2" 

9 2 g 2 

We see perturbatively in g that x(%) ^ s a massless field. 



e 9X - l - g x - - g 2 x 2 



(57) 



Concerning the boundary conditions, we must consider the system inside a large sphere of radius R (10 4 
10 -2 pc < R< 100 pc). That is, all integrals are computed over such large sphere. 



Using cq.(|30|) the particle density takes now the form 



p(r) 



1 



(0 = -; 



V 2 x(r*) 



H 2 5 r 



''/J ''/J T eff 

It is convenient to renormalize the particle density by its stationary value 6 — 



i - gx(r) 



P{r)ren = J p{r) = 

5 g 2 



(58) 



We see that in the S — > limit the interaction in eq.(p7|) vanishes. No infrared divergences appear in the Feynman 
graphs calculations, since we work on a very large but finite volume of size R. Hence, in the S — ► limit, the whole 
perturbation series around 4>o reduces to the zeroth order term. 
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The constant saddle point <po fails to catch up the whole field theory content. In fact, more information arises 
perturbing around the stationary point (f> c (r) given by eq.(^) 

Using eqs.(Bq), ([43]), (p7|) and (pq) we obtain for the density correlator in the 5^0 limit, 



2D 



C(fi,r 2 ) = —r { exp 



r 



9 



C D {p \n - f 2 \) 



D-2 



C D (fi |fi - r 2 |) 



D-2 



For large distances, we find 



\rt— F 2 |— »°o 



2Cf, In-fal 2 ^- 2 ) 



0(|ri-r 2 | 



-3(D-2) 



(59) 



That is, the 0-field theory scales. Namely, the theory behaves critically for a continuum set of values of fj, and 

T eff- 

Notice that the density correlator C(fi,r2) behaves for large distances as the correlator of x(^) 2 - This stems from 
the fact that x(r) 2 is the most relevant operator in the series expansion of the density (pq) 



P(r)r 



1 



M x(rY + 0( X 6 ) 



(60) 



As remarked above, the constant stationary point 4>q = log S — ► — oo only produces the zeroth order of perturbation 
theory. More information arises perturbing around the stationary point 4> c {r) given by eq.(|5^) 63). 



C. Renormalization Group Finite Size Scaling Analysis 



As is well known |12|-|14|, physical quantities for infinite volume systems diverge at the critical point as A to a 
negative power. A measures the distance to the criti cal poi nt. (In condensed matter and spin systems, A is proportional 
to the temperature minus the critical temperature |l^,|lj]). One has for the correlation length £, 



and for the specific heat (per unit volume) C, 



£(A) ~ A"" 



C - A- a 



(61) 

(which in spin systems describes the magne- 



Correlation functions scale at criticality. For example, the scalar field 
tization) scales as, 

< 4>(r)cj>{0) >~ r~ l ~ v . 

The critical exponents v, a and r] are pure numbers that depend only on the universality class |l2|-[l4| . 

For a finite volume system, all physical quantities are finite at the critical point. Indeed, for a system whose size 
R is large, the physical magnitudes take large values at the critical point. Thus, for large R, one can use the infinite 
volume theory to treat finite size systems at criticality. In particular, the correlation length provides the relevant 
physical length £ ~ R. This implies that 



A - R 



-l/v 



(62) 



We can apply these concepts to the (/)-theory since, as we have seen in the previous section, it exhibits scaling in a finite 
volume ~ R 3 . Namely, the two points correlation function exhibits a power-like behaviour in perturbation theory as 
shown by eq.([59"|). This happens for a continuum set of values of T e // arid /x 2 . Therefore, changing /j, 2 /T e ff keeps 
the theory in the scaling region. At the point /j 2 /T e ff = 0, the partition function Z is singular. From eq.d22), we 
shall thus identify 



A = 



3/2 



T, 



eff 



V 2vr / 



(63) 



Notice that the critical point A = 0, corresponds to zero fugacity. 
Thus, the partition function in the scaling regime can be written as 
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Z(A) = 



V<f> e 



-S'+A J d D x e* (£ > 



where S* stands for the action (|23|) at the critical point A = 0. 
We define the renormalized mass density as 

mp(x) ren = m e^^' 



(64) 



(65) 



and we identify it with the energy density in the renormalization group. [Also called the 'thermal perturbation 
operator']. This identification follows from the fact that they are the most relevant positive definite operators. 
Moreover, such identification is supported by the perturbative result (|60|). 
In the scaling regime we have ]13| for the logarithm of the partition function 



1 K 

- logZ(A) = 

V (2 — a)[l — a) 

where F(A) is an analytic function of A around the origin 



A 



F(A) 



(66) 



F(A) =F + aA+-bA 2 + ... . 

V = R D stands for the volume and Fq, K, a and b are constants. 

Calculating the logarithmic derivative of Z(A) with respect to A from eqs. 
results yields 



lAlogZ(A)= a +-^A 1 -« 
V OA 1 — a 



vl 



d x < e 



4,(3) 



and from (66) and equating the 
> • (67) 



where we used the scaling relation a = 2 — vD [p"3 14 . 

We can apply here finite size scaling arguments and replace A by ~ R~~ [eq. (|62|)] , 

^-io g z(A) = y a + -^i? 1 /- + ... . 

a A 1 — a 

Recalling eq. ([35]) , we can express the mass contained in a region of size R as 



M(R) = m 



3 *<*> d D x 



(68) 



Using eq.(|67|) we find 



and 



< M(R) >=mVa + m 



K 



1 — a 



R- 



< p(r) >= ma + m 



K 



■ v D for r of order ~ R. 



where ilo is the surface of the unit sphere in Z?-dimensions. 

The energy density correlation function is known in general in the scaling region (see refs. |l3| 
therefore write for the density-density correlators (|34j) in D space dimensions 



C(fi,r 2 ) ~ \n - r 2 \ 



--2D 



(69) 

|). We can 
(70) 



where both fi and f*2 are inside the finite volume ~ R D . 

The perturbative calculation ( ^9[ ) matches with this result for v — i. That is, the mean field value for the exponent 

v. 

Let us now compute the second derivative of log-Z(A) with respect to A in two ways. We find from eq.(|6l 

^logZ(A) = V[A- a K + b+...] . 
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We get from eq.(p4|), 



(71) 



where we used eq.(|62|), eq.(|70|) and the scaling relation a = 2 — OJTJ]. We conclude that the scaling behaviours, 
eq.(|6^) for the partition function, eq. (pT|) for the specific heat and eq. (|70|) for the two points correlator are consistent. 
In addition, eqs.(p8|) and (|7l|) yield for the mass fluctuations squared 



(AM(R)) 2 = < M 2 > — < M > 2 ~ J d D x d D y C(x,y) ~ R 2d " . 

Hence, 

AM(R) ~ R dH . (72) 

The scaling exponent v can be identified with the inverse Haussdorf (fractal) dimension dn of the system 

d H = - ■ 
v 

In this way, AM ~ R da according to the usual definition of fractal dimensions . 

Using eq. (|7p| ) we can compute the average potential energy in three space dimensions as 

<V>=\pGm 2 I" C(x,ff) ~ Rl- 1 . 



2 J\S-y\>a \x-y\ 



From here and eq.(72) we get as virial estimate for the atoms kinetic energy 

< V > 
< AM(R) > 



, < V > i , 

< V 2 >= —7^ R-^ 



This corresponds to a velocity dispersion 

Au~ J R3(i- 1 ) . (73) 
That is, we predict [see eq.p)] a new scaling relation 



2\v J 2 

The calculation of the critical amplitudes [that is, the coefficients in front of the powers of R in eqs.(|70|), ( |72|) and 
(f73|)] is beyond the scope of the present paper p^ ]. 

D. Values of the scaling exponents and the fractal dimensions 

The scaling exponents is, a considered in sec IIIC can be computed through the renormalization group approach. 
The case of a single component scalar field has been extensively studied in the literature |T^,[l4|,[l7j . Very probably, 
there is an unique, infrared stable fixed point in three space dimensions: the Ising model fixed point. Such non- 
perturbative fixed point is reached in the long scale regime independently of the initial shape of the interaction u(4>) 
[eq.©] 0. 

The numerical values of the scaling exponents associated to the Ising model fixed point are 

v = 0.631... , d H = 1.585 .. . , r) = 0.037... and a = 0.107... . (74) 

In the cj> field model there are two dimcnsionful parameters: /i and T e ff. The dimensionless combination 
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^3/2 ™15/4 

<? = pT eff = (8*)"* V* T ™ (75) 

acts as the coupling constant for the non-linear fluctuations of the held <j>. 

Let us consider a gas formed by neutral hydrogen at thermal equilibrium with the cosmic microwave background. 
We set T = 2.73 K and estimate the fugacity z using the ideal gas value 



(-) 



3/2 

P 



Here we use p = So atoms cm 3 for the ISM density and So — 10 10 . Eq. eqs.(|22|) yields 
1 1 



2.7 -^= AU - 30 AU and g 2 = [i T ef f = 4.9 10" 58 y/So ~ 5 10~ 53 . (76) 

This extremely low value for g 2 suggests that the perturbative calculation [sec. IIIB] may apply here yielding the 
mean held values for the exponents, i. e. 

v = 1/2 , d H = 2 , 7] = and a = . (77) 

That is, the effective coupling constant grows with the scale according to the renormalization group how (towards 
the Ising fixed point). Now, if the extremely low value of the initial coupling eq.([76[) applies, the perturbative result 
(mean held) will hold for many scales (the effective g grows roughly as the length). 

/i -1 indicates the order of the smallest distance where the scaling regime applies. A safe lower bound supported by 
observations is around 20 AU ~ 3. 10 14 cm , in agreement with our estimate. 

Our theoretical predictions for AM(R) and Av [eqs.([72|) and (f73j)] both for the Ising cq.(f74|) and for the mean held 
values eq.(Q), are in agreement with the astronomical observations [eq. (jl|)] - The present observational bounds on the 
data are larger than the difference between the mean held and Ising values of the exponents dn and q. 

Further theoretical work in the 0-theory will determine whether the scaling behaviour is given by the mean held or 
by the Ising hxed point |p2fl . 



E. The two dimensional gas and random surfaces fractal dimensions 



In the two dimensional case {D — 2) the partition function (|45|) describes the Liouville model that arises in string 
theory [ pp[ and in the theory of random surfaces (also called two-dimensional quantum gravity). For strings in 
c-dimensional Euclidean space the partition function takes the form pOfl 

Z e = J J V4> e"^ 2 / d2 4iW) 2 +m 2 e««] _ (7g) 

This coincides with eq.(|45|) at D — 2 provided we hip the sign of /j, 2 and identify the parameters (^) as follows, 

T = Gm 2^_£ ; ^ = z Gm 3 . (79) 

Ref. jn| states that dn = 4 for c < 1, dn — 3 for c = 2 and dn = 2 for c > 4. In our context this means 

25 11 
d H = 2 for T < — Gm 2 , d H = 3 for T = 2 Gm 2 and d H = 4 for T > — Gm 2 . 

For c — > oo, g 2 — > and we can use the perturbative result ( |59| ) yielding v = i dn = 2 in agreement with the 
above discussion for c > 4. 
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F. Stationary points and the Jeans length 

The stationary points of the 0-field partition function ( |2l| ) are given by the non-linear partial differential equation 

V 2 </> = —p 2 e*<® . 

In terms of the gravitational potential U(x) [see eq. (p3[)], this takes the form 



V 2 C/(r) = AttGz'i 



( mT 



3/2 



,-fc/(fO 



(80) 



This corresponds to the Poisson equation for a thermal matter distribution fulfilling an ideal gas in hydrostatic 
equilibrium, as can be seen as follows |l6| . The hydrostatic equilibrium condition 

VP(f) = -mp(f) VU(r) , 

where P(r) stands for the pressure, combined with the equation of state for the ideal gas 

P = Tp, 

yields for the particle density 

p(r)=p e-^ u ^ , 
where pa is a constant. Inserting this relation into the Poisson equation 

V 2 [/(r) = 4wGmp{r) 

yields eq. (|80|) with 



po 



mT 
~2T 



3/2 



For large r, eq.(80) gives a density decaying as 



p{r) 



T 1 



2nGm r 2 



U{r) r =°° — log 
m 



2ttGpq 
T 



(81) 



(82) 



Notice that this density, which describes a single stationary solution, decays for large r faster than the density (|6_ 
governed by thermal fluctuations. 

Spherically symmetric solutions of eq.([30|) has been studied in detail |18|j . The small fluctuations around such 
isothermal spherical solutions as well as the stability problem were studied in |19| . 



The Jeans distance is in this context, 



m ytGrnpo ~ VlG m 7 / 4 T 1 / 4 



(83) 



This distance precisely coincides with pT 1 [see eq.(^2])] up to an inessential numerical coefficient (y/12/n). Hence, p, 
the only dimensionful parameter in the 0-theory can be interpreted as the inverse of the Jeans distance. 
We want to notice that in the critical regime, dj grows as 



R 



d H /2 



(84) 



since po = A ~ Rr da vanishes as can be seen from eqs.(62), ( |6^ ) and (|83|). In this tree level estimate we should use 
for consistency the mean field value dn — 2, which yields dj ~ R. 

This shows that the Jeans distance is of the order of the size of the system. The Jeans distance scales and the 
instability is therefore present for all sizes R. 

Had dj being of order larger than R, the Jeans instability would be absent. 



1G 



The fact that the Jeans instability is present precisely at dj ~ R is probably essential to the scaling regime and 
to the self-similar (fractal) structure of the gravitational gas. 

The dimensionless coupling constant g 2 can be written from eqs.(|63|) and (75) as 



2m 



ttG 



A 



Hence, the tree level coupling scales as 
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R- 1 



Direct perturbative calculations explicitly exhibit such scaling behaviour 
We can express g 2 in terms of dj and po as follows, 



(12tt) 3 / 2 
Po d 3 j 



ir 2 p 3 
Po 



This shows that g 2 is, at the tree level, the inverse of the number of particles inside a Jeans volume. 

Eq.(^) applies to the tree level Jeans length or tree level pT 1 ■ We can furthermore estimate the Jeans length using 
the renormalization group behaviour of the physical quantities derived in sec. III.C. Setting, 



< dj >= 



< Av > 



y/Gm < Ap> ' 



we find from eqs.(|69[) and ([j 



< dj >~ R 



Namely, we find again that the Jeans length grows as the size R. 



IV. DISCUSSION 



In previous sections we ignored gravitational forces external to the gas like stars etc. Adding a fixed external mass 
density p e xt(r) amounts to introduce an external source 

J(r) = -T eff p ext (r) , 

in cq.(j2^). Such term will obviously affect correlation functions, the mass density, etc. except when we look at the 
scaling behaviour which is governed by the critical point. That is, the values we find for the scaling exponents du 
and q are stable under external perturbations. 

We considered all atoms with the same mass in the gravitational gas. It is easy to generalize the transformation 
into the 0-field presented in section II for a mixture of several kinds of atoms. Let us consider n species of atoms with 
masses m a , 1 < a < n. Repeating the steps from eq.(|T^) to (^j|) yields again a field theory with a single scalar field 
but the action now takes the form 



l 



T, 



<ff 



d 3 3 



Ma e " 



■Hi 



(85) 



where 



z a G m 3 / 2 m 2 Vf 



and m is just a reference mass. 

Correlation functions, mass densities and other observables will obviously depend on the number of species, their 
masses and fugacities but it is easy to see that the fixed points and scaling exponents are exactly the same as for the 
0-field theory ©-(§!). 
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We want to notice that there is an important difference between the behaviour of the gravitational gas and the 
spin models (and all other statistical models in the same universality class). For the gravitational gas we find scaling 
behaviour for a full range of temperatures and couplings. For spin models scaling only appears at the critical value 
of the temperature. At the critical temperature the correlation length £ is infinite and the theory is massless. For 
temperatures near the critical one, i. e. in the critical domain, £ is finite (although very large compared with the 
lattice spacing) and the correlation functions decrease as ~ e~ r ^ for large distances r. Fluctuations of the relevant 
operators support perturbations which can be interpreted as massive excitations. Such (massive) behaviour does not 
appear for the gravitational gas. The ISM correlators scale exhibiting power-law behaviour. This feature is connected 
with the scale invariant character of the Newtonian force and its infinite range. 

The hypothesis of strict thermal equilibrium does not apply to the ISM as a whole where temperatures range from 
5 to 50 K and even 1000 K. However, since the scaling behaviour is independent of the temperature, it applies to 
each region of the ISM in thermal equilibrium. Therefore, our theory applies provided thermal equilibrium holds in 
regions or clouds. 

We have developed here the theory of a gravitationally interacting ensemble of bodies at a fixed temperature. In 
a real situation like the ISM, gravitational perturbations from external masses, as well as other perturbations are 
present. We have shown that the scaling solution is stable with respect to the gravitational perturbations. It is well 
known that solutions based on a fixed point are generally quite robust. 

Our theory especially applies to the interstellar medium far from star forming regions, which can be locally far 
from thermal equilibrium, and where ionised gas at 10 4 K together with coronal gas at 10 6 K can coexist with the cold 
interstellar medium. In the outer parts of galaxies, devoid of star formation, the ideal isothermal conditions arc met 
0. Inside the Galaxy, large regions satisfy also the near isothermal criterium, and these are precisely the regions 
where scaling laws are the best verified. Globally over the Galaxy, the fraction of the gas in the hot ionised phase 
represents a negligible mass, a few percents, although occupying a significant volume. Hence, this hot ionised gas is 
a perturbation which may not change the fixed point behaviour of the thermal self-gravitating gas. 

In ref. pij a connection between a gravitational gas of galaxies in an expanding universe and the Ising model is 
conjectured. However, the unproven identification made in ref. p4| of the mass density contrast with the Ising spin 
leads to scaling exponents different from ours. 
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